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Abstract. A model of the grain surface chemistry involving the accretion of atoms of two different elements, X and 
Y, and their reactions to form species X2, XY, and Y2 was examined for a wide range of choices for the values of its 
£f~) ' three free parameters - the accretion rate of X and Y, the desorption rate of X and the grain surface sweeping time 

of Y, all considered relative to the grain surface sweeping rate of X. Relative production rates of the diatomics were 
calculated with five methods involving, respectively, a high-order truncation of the master equation, a low-order 
truncation of the master equation, the standard deterministic rate equation approach, a modified rate equation 
approach and a set of approximations which are in some cases appropriate for accretion dominated chemistry. The 
accuracies of the relative production rates calculated with the different methods were assessed for the wide range 
of model parameters. The more accurate of the low-truncation master equation calculations and the standard 
deterministic rate equation approach gives results which are in most cases within ten or twenty per cent of the 
' results given by the high-truncation master equation calculations. For many cases, the more accurate of the low 

order truncation and the standard deterministic rate equation approaches is indicated by a consideration of the 
■ average number of atoms of the two species on the grain's surface. 

^p^i Key words, astrochemistry - molecular processes - ISM: clouds - dust, extinction - ISM: molecules 
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^ ; 1. Introduction Robinson (1977) and Tielens & Hagen (1982). The model 

of Allen & Robinson (1977) was constructed under the 
> . If species accrete from the gas phase onto the surfaces of assumption that when a mo i eC ule is produced in a grain 
; grains more quickly than they diffuse across a grain surface surface reactioil] it i mme diately enters the gas phase; the 

' and react with each 0ther ' the Surface chemistr y is said t0 model was extended by Tielens & Hagen (1982) to al- 
C3 ; be taking place in the reaction limit. Such chemistry has bw for ^ accumulation of grain mantles . In their mode l, 
been treated by Pickles & Williams (1977), who used a Tieleng & Hagen (1Q82) used Monte Carb techniques 
system of rate equations similar to those employed in gas- tQ calculate the s t ea dy-state concentrations of mantled 
phase chemistry. This approach is often referred to as the specieg _ Laterj yarious Qther authorg constructed Monte 
standard deterministic rate equation method. Carlo q{ gm[n surface chemistry (e _ g _ Tielens & 
In some cases, species diffuse across grain surfaces and Allamandola 1987; Tielens 1995 (unpublished); Charnley 
react with each other faster than they accrete. Therefore, et al 1997; T i e lens & Charnley 1997; Charney 2001). 
when a species is accreted, and there is already another charnley (1998) used Monte Carlo methods to solve the 
reactive species on the grain, they are likely to react before master equation (which gives the probability that there 
another particle accretes. The rates of reaction are there- j s a cer t, am number of particles of each species) govern- 
fore limited by the rates of accretion; the chemistry is said ing gas _ phase chemistry, and suggested that a similar ap- 
to be occurring in the accretion limit. Under these condi- proach may be ap pii ca ble to grain surface chemistry, 
tions, the average numbers of reactive species on a grain 
are small, and a stochastic model is required. The first 

attempts at stochastic modelling of interstellar grain sur- As Monte Carl ° simulations are computationally too 

face chemistry in the accretion limit were made by Allen & expensive to use in the study of large systems, Caselli et 

al. (1998; 2002) and Shalabiea et al. (1998) introduced 

Send offprint requests to: J.G.L. Rae e-mail: semi-empirical modifications to the standard determinin- 

jglr@ast.leeds.ac.uk stic rate equations, with the aim of solving the problem 



2 



J.G.L. Rae et al.: Diffusive grain surface chemistry 



of inaccuracy of rate equation methods in the accretion 
limit. 

Biham et al. (2001) and Green et al. (2001) studied, 
with a master equation approach, the production rate of 
H 2 on surfaces. Green et al. (2001) used a generating 
function method to obtain an analytic solution for the 
H2 production rate. They also investigated more complex 
chemistries, using sparse matrix techniques to obtain the 
probabilities that a grain contains given numbers of par- 
ticles of each species. With this approach, the calculation 
must be truncated at a certain maximum number of parti- 
cles of each species, i.e. there is some truncation value for 
the number of particles, above which the probabilities are 
assumed to be zero. With a high enough truncation value, 
the results are almost exact. Green et al. (2001) used the 
probabilities obtained in this way to calculate the rates 
of production of molecules on grain surfaces. They solved 
the master equation numerically for the system of H and 
O, reacting on grain surfaces to form H2, O2, and OH, 
and also for the system of H, O, and N, reacting to form 
H 2 , 2 , OH, NH, N 2 , NO, N0 2 , H 2 0, and NH 3 . They 
compared the results given by the solution of the master 
equation with a high truncation value, to those given by 
(i) the use of the master equation method with a low trun- 
cation value; (ii) the standard deterministic rate equation; 
and (iii) an approximate method, based on the assump- 
tion that the sweeping rate of atomic hydrogen on grain 
surfaces is so great that a reactive species already on the 
surface of a grain will react as soon as an H atom is ac- 
creted. 

Stantcheva et al. (2002) solved the master equation for 
the system of H, O, and CO reacting on the surfaces of 
grains to produce 2 , H 2 , H 2 0, C0 2 , H 2 CO, and CH 3 OH. 
They investigated the accuracy of results given by the 
method for different truncation values and also the ac- 
curacies of results given by other approaches. 

In this paper, we examine the system of atoms X and 
Y reacting to form X2, XY, and Y2 on grain surfaces. 
We study a wide variety of adsorption, desorption, and 
diffusion rates, covering the accretion limit, the reaction 
limit, and cases intermediate between them. We compare 
the almost exact results, obtained through the solution 
of the master equation for high truncation values, with 
the results obtained (i) through the solution of the master 
equation with low truncation values; (ii) with the approx- 
imate method used by Green et al. (2001); (iii) with the 
standard deterministic rate equation approach; and (iv) 
with a modified rate equation approach similar to that 
employed by Caselli et al. (2002). Although the network 
we study is smaller than that of Stantcheva et al. (2002), 
we examine a much larger range of adsorption, desorp- 
tion, and diffusion rates, thereby gaining insight into the 
parameter ranges in which each method may be applied. 



2. The various approaches to the problem 

We consider a system in which atoms X and Y accrete 
from the gas phase onto the surfaces of grains, diffuse 
across the grain surfaces, and react with each other to 
form X2, Y2, and XY, i.e. in which the following reactions 
take place: 



x + x -> x 2 

Y + X -> XY 

Y + Y -> Y 2 . 



(1) 
(2) 
(3) 



We assume that the X 2 , XY, and Y2 molecules remain on 
the grain surfaces, but that they do not react with each 
other or with X or Y. We allow for desorption of X from 
grains, but assume that desorption of Y is negligible. Let 
ax and ay be the rates of accretion of X and Y respec- 
tively, and dx be the rate of desorption of X. Let f s (X) 
and t s (Y) be the grain surface sweeping times for X and 
Y respectively. To reduce the number of free parameters, 
we assume 



ax = ay 



(4) 



This assumption is likely to be reasonable at some den- 
sities in dense cores if X is hydrogen and Y is oxy- 
gen. Furthermore, in a number of previous studies (e.g. 
Stantcheva et al. 2002), the ratio of accretion rates has 
been varied. 

We normalise with respect to t s (X) as follows. 



a = a x i s (X) 
6 = d x t s {X) 

= uoo 

t s (X) 

The dimensionless rates for Reactions (Q to © are 

kxx = 2 

kxy = 1 + t^ 1 

Kyy = 2t~ 1 



(5) 
(6) 

(7) 



(8) 
(9) 
(10) 



Eqs. iJSJ to HI L)| l may be compared to the rates given by 
Eq. (9) of Caselli et al. (1998) for two arbitrary species I 
and J. 

The production rates of X2, XY, and Y2 were calcu- 
lated with the various approaches for 13 x 13 x 13 sets 
of values of a,S,r between 0.1 and 10. Calculations were 
also performed for 9 x 13 x 9 sets of 0.01 < a < 0.09, 
0.1 < S < 10.0, 0.01 < t < 0.09. 

We now discuss in detail the master equation ap- 
proach, the standard and modified rate equation ap- 
proaches, and the approximate method of Green et al. 
(2001). 

2.1. Master equation approach 



In Sect. |2 the various approaches to the problem are Green et al. (2001) gave the master equation for the hy- 
discussed. We present our results in Sect. and Sect. 01 drogen and oxygen system in their Eq. (29), and the rates 
concludes the paper. of production of H2, OH, and O2 in their Eqs. (26) to (28). 
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The generalisation of their equations to species X, Y, X 2 , 
XY, and Y 2 is trivial. The average numbers of X and Y 
atoms on the surface of a grain are: 



(N(X)) = jTiP(i,j) 



i=0 

3=0 



(N(Y)) = X)jP(i,j) 



(11) 
(12) 



«=0 

3=0 



where P(i,j) is the probability that on the surface of a 
grain there are i atoms of species X and j atoms of species 
Y. We assumed the system to be in steady state, so that 
dP ^ J - ) = for all and followed Green et al. (2001) 
in solving the master equation through the inversion of a 
sparse matrix by the row-indexed method. As discussed 
in Scct.^ with this approach it is necessesary to truncate 
the sums at particular values of i and j. They were trun- 
cated first at i — j — 2 (which will be referred to as the 
low truncation case) , and then at higher values (which will 
be referred to as the high truncation case), above which 
the results did not change appreciably if i and j were in- 
creased further. The results given by the master equation 
approach with high truncation values of i and j - typi- 
cally i = j > 5 - can be considered to be exact, and will 
be referred to as the exact results. 



2.2. Standard rate equation approach 

The standard deterministic rate equations for the H, O, 
H2, OH, and O2 system are given by Caselli et al. (1998) 
in their Eqs. (4) to (8). They can be generalised easily to 
the X, Y, X2, XY, and Y2 system. The rate equations for 
X and Y were integrated with a Gear algorithm until the 
system reached steady state. Then, the production rates of 
X2, XY, and Y2 were calculated from their rate equations. 



2.3. Modified rate equation approach 

We introduced modifications to the rate equations, similar 
to those suggested by Caselli et al. (1998, 2002) in their 
attempt to develop a set of deterministic equations appro- 
priate in both the reaction limit and the accretion limit. 
The modifications we used were based on those used by 
Caselli et al. (2002). 

We define Q\ and 9 2 as 



72 = 



K X ,x(iV(X)) + «x,Y(iV(Y)) 
a 

Ky , y (N{Y))+k x , y (N{X)) 



a 



(13) 
(14) 



In cases in which 6\ < 1 and 9 2 < 1, the standard deter- 
ministic rate equations were used. Otherwise, the equa- 
tions were modified as follows. 

In the case that 9\ > 1 and 9 2 < 1, the probability 
that species Y reacts is either less than or equal to 1, 
so that for Y2 formation the standard deterministic rate 



equation approach can be used. However the probability 
that species X reacts is greater than 1, and the equations 
for the production of X2 and XY must be modified. The 
rate coefficients k x ,x and kx,y are replaced by the larger 
of a and S (Caselli et al. 1998, 2002); this quantity will 
be denoted by k' x . In addition, the formation rate of X2 is 
multiplied by the probability that X reacts with another 
X instead of with a Y, and the rate of formation of XY 
is multiplied by the probability that X reacts with Y and 
not with another X (following Caselli et al. 2002). The 
production rates, T(X 2 ) and T(XY), of X 2 and XY, in 
this approach become 



r(x 2 ) = k' x (n(x)) 
r(XY) = k' x (n(x)) 



k x ^(N(X)) 



«x,x(^(X))+kx,y<^(Y)) 

K X ,Y(iV(Y)) 
kx,y<JV(Y)>+k XiX (JV(X)> 



(15) 
(16) 



After l|15|) , Ijltjfl , and the standard rate equation for Y 2 are 
used, 02 is recalculated; if it is found to exceed unity then 
the approach described below for cases in which 9\ > 1 
and 9 2 > 1 is applied. 

The case in which 9\ < 1 and 9 2 > 1 is the same as 
the previous case, except that X is replaced by Y and vice 
versa. So T(X 2 ) can be calculated with the standard deter- 
mistic rate equation approach, and the rates of formation, 
r(XY) and T(Y 2 ), of XY and Y 2 , are 



r(XY) = k' y (n(y)) 

T(Y 2 ) = k y {N(Y)) 



k x ,y(N(X)) 



KX,Y(AT(X))+K Y ,Y<iV(Y)) 
K Y ,Y(iV(Y)) 

k y ,y(N(Y))+ Kx .y(N(X)) 



(17) 
(18) 



where Kyis the larger of the accretion rate and the des- 
orption rate of Y; however because we are neglecting des- 
orption of Y, k' y = a. After Eqs. (fTTjl and (fT£)) and the 
standard rate equation for X 2 production are used, 9\ is 
recalculated; if it is found to be greater than unity then 
the approach described below for cases in which 9\ > 1 
and 92 > 1 is applied. 

In the case in which both 9\ and 92 are greater than 1, 
the rate equations for X2, XY, and Y 2 must be modified. 
For the formation rates of X 2 and Y 2 we used Eqs. <|15fl 
and (|18f) . respectively. For the XY formation rate, we used 
Eq. (HI) when 9\ > 6> 2 , and Eq. (JTJl when 0\ < 9 2 . 

In all cases, the rates of change of (N(X)) and (N(Y)) 

are 



d(iV(X)) 
dt 

d(N(Y)) 
dt 



a - 5{N(X)) - 2T(X 2 ) - T(XY) (19) 
a - 2r(Y 2 ) - T(XY) (20) 



Eqs. (ftT))) and f2T)| were integrated with a Gear algo- 
rithm until the system reached steady state. At each step 
in the integration, the rates of production of X2 , XY, and 
Y2 were calculated using the standard deterministic, or 
modified, rate equations as appropriate, and were used in 



the calculations of 



*(N(X)) 



dt 



and 



*(N(Y)) 

dt 
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2.4. Approximate Method 

Green et al. (2001) introduced an analytic approximation 
to the master equation method, based on the assumption 
that the sweeping rate of X is so fast that if there is a 
reactive species (X or Y in this case) on the surface of a 
grain, a X atom accreting onto the grain will react im- 
mediately. This is equivalent to assuming that kx,x 3> a, 
and that the probability of there being more than one re- 
active species on the surface of a grain is small. Green et 
al. (2001) showed that, under these assumptions 



r(x 2 ) 
r(XY) 
r(Y 2 ) 

where 
P(0,0) 



a 2 P(0,0) 



2a 

,2 



8 

1 | 1 
2a + 5 2a 
Ky y ol 



ky y + 2a 2 



P(0,0) 
P(0,0) 



2a + S 



(21) 
(22) 
(23) 

(24) 



3. Results and discussion 

For each of the methods described in Sect. [3 the produc- 
tion rates of X 2 , XY, and Y 2 were evaluated with a, 8, 
kx,x, «x,y, and ky,y given by Eqs. J5j to l|10|) . with a, <5, 
and t having various values, as described in Sect.HJ 

Following Caselli et al. (1998), we calculated the rel- 
ative abundances of X 2 , XY, and Y 2 on the surfaces of 
grains as 



x(H 2 ) 
x(XY) 
x(0 2 ) 



r(x 2 ) 


r(x 2 ) 


4- r(XY) - 


^r(Y 2 ) 




r(XY) 




r(x 2 ) 


4- T(XY) - 


i- r(Y 2 ) 




r(Y 2 ) 




r(x 2 ) 


4- r(XY) - 


h r(Y 2 ) 



(25) 
(26) 
(27) 



We calculated the difference between the exact results, 
and the results obtained (i) through the solution of the 
master equation for the low-truncation case; (ii) with the 
standard rate equation approach; (hi) with the modified 
rate equation approach; and (iv) with the approximate 
method. The results are plotted against a and 8, for r = 
1.0, in Figs. [J to H 

Figs, n to 01 show the magnitudes of the percentage 
discrepancies between: the results given by the standard 
rate equation approach, and the exact results (Fig. 
the results given by the modified rate equation approach, 
and the exact results (Fig. [2); the results given by the 
approximate method, and the exact results (Fig. EJ; and 
the results given by the master equation approach for the 
low truncation case, and the exact results (Fig.EJ- 

Due to lack of space, the results are plotted against a 
and 8 for r = 1.0 only. Results for other values of r can 
be found at http : //ast . leeds . ac . uk/^jglr. In each of 



Figs.[I]to|H plot (a) shows the discrepancies in x(X 2 ); plot 

(b) shows the discrepancies in x(XY); and plot (c) shows 
the discrepancies in x(Y 2 ). 

In Fig. the percentage discrepancies are plotted 
against (N(X)) and (N(Y)), the average numbers of X and 
Y atoms on the surface of a grain, as given by the exact 
method and Eqs. I|ll(l and (|12f> . The white gaps at the top 
left-hand and bottom left-hand corners of the plots in Fig. 
El are caused by the fact that we did not explore regions 
of a-S-T space that yielded these values of (JV(X)) and 
(N(Y)). The uneven contours in Fig. are caused by the 
fact that any point in (N(X))-{N(Y)) space does not cor- 
respond to a single point in a-S-T space, and therefore the 
discrepancies do not depend only on (iV(X)) and (N(Y)). 
In Fig.Elthe discrepancies are plotted against (N(X)) and 
(N(Y)) for r = 1.0 only. Again the white areas are regions 
of (N(X))-(N(Y)) space to which no point in a-8 space 
corresponds for r = 1.0. In Figs. El and El plots (a), (b) 
and (c) are for the standard rate equation approach, plots 
(d), (e), and (f) for the modified rate equation approach, 
plots (g), (h), and (i) for the approximate method, and 
plots (j), (k), and (1) for the low-truncation case of the 
master equation approach. Plots (a), (d), (g), and (j) give 
the percentage discrepancy in x(X 2 ), plots (b), (e), (h), 
and (k) the percentage discrepancy in x(XY), and plots 

(c) , (f), (i), and (1) the percentage discrepancy in x(Y 2 ). 
In Figs. Eland© the ranges of (N (X))-(N (Y)) are differ- 
ent because there is no desorption of Y, so X and Y are 
not symmetric. 

In all of Figs. ^ to El lighter regions indicate lower 
percentage discrepancies, and darker regions higher per- 
centage discrepancies. 

Although only the results for r = 1.0 are shown in Fig. 
^ it was found that the results given by the rate equation 
approach are inaccurate for the case of low a, high 8, and 
low r, corresponding to situations in which the rates of 
accretion are low and the rates of X desorption, and of 
Reactions J2J and J3J, are high, so that the average surface 
population on a grain is low. This conslusion is confirmed 
by Figs.EK _c andE^-c, which show that this approach is 
least accurate for small (N(X)) and (N(Y)). These are the 
conditions under which it is known that the rate equation 
approach breaks down. 

The inaccuracies in the results given by the modified 
rate equation approach are shown in Fig. I^for r = 1.0, 
and in Figs.Eli-f andEtf-f- In most cases these results are 
no more accurate than those of the standard rate equa- 
tion approach, and sometimes the accuracy is much less. 
Stantcheva et al. (2002) studied the system of H and O 
reacting on the surfaces of grains to form H 2 , OH, and 
2 , and gave in their Figs. 1 to 3 results for the values 
of x(H 2 ), x(OH), and x(0 2 ) calculated by various dif- 
ferent methods including the modified rate equation ap- 
proach and the master equation method. The modified 
rate equation approach they used was the same as that 
used by Stantcheva et al. (2001), which was based on that 
used by Caselli et al. (1998). They found good agreement 
between results obtained with these two approaches for 




Fig. 1. Percentage discrepancies in results given by stan- 
dard rate equation approach for r = 1.0. a: discrepancy in 
x(X2); b: discrepancy in s(XY); c: discrepancy in x(Y2). 
Darker regions indicate greater discrepancies. 



Fig. 2. Percentage discrepancies in results given by mod- 
ified rate equation approach for r = 1.0. a: discrepancy in 
x(X.2); b: discrepancy in i(XY); c: discrepancy in x(Y2). 
Darker regions indicate greater discrepancies. 



the parameters they considered, which were mostly rather ucs for ax, ay, dx, t s (X), and i s (Y)as Hasegawa et al. 
different from those used by us. They used the same val- (1992), Caselli et al. (1998), and Green et al. (2001). The 
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Fig. 3. Percentage discrepancies in results given by ap- Fig. 4. Percentage discrepancies in results given by master 

proximate method for r = 1.0. a: discrepancy in x(X.2); equation approach with low truncation values for r = 1.0. 

b: discrepancy in x(XY); c: discrepancy in a;(Y 2 ). Darker a: discrepancy in x(X 2 ); b: discrepancy in i(XY); c: dis- 

regions indicate greater discrepancies. crepancy in a;(Y 2 ). Darker regions indicate greater dis- 
crepancies. 

accretion rate they used for oxygen was different from 
that used for hydrogen. Their parameters were equiv- 
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alent to 1CT 12 < a H i s (H) < 1(T 5 ; a £ s (H) - 10" 10 ; 
cfe* s ( H ) - 1CT 8 ; and M|| - 10 9 . This corresponds to 
a region below and to the left of that shown in the plots 
in Fig. [3 and for a higher value of r. It can be seen that 
towards the bottom left-hand corner of Fig. EH, the results 
are more accurate than elsewhere, which is consistent with 
the findings of Stantcheva et al. (2002). 

Fig. 13 shows that the approximate method gives inac- 
curate results for x(X.2) when 8 is high. If the desorption 
rate is high, the assumption that an X atom will always 
react if it accretes onto a grain that is already populated 
is invalid. It was also found that there are inaccuracies in 
the values calculated for ir(Y 2 ) for high a, for all values 
of t. In these cases, the accretion rate is comparable to, 
or greater than, the rate of reaction of Y, and the popula- 
tion of grains containing several Y atoms is non-negligible. 
At high values of r, the results were found to be inaccu- 
rate for a wide range of values of a and 5, particularly for 
x(X.2) and a;(Y 2 ); x(XY) is only inaccurate for high 6 and 
low a. The approximate equations were constructed un- 
der the assumption that if there are two atoms of species 
Y on the surface of a grain, and a third accretes, Y 2 will 
not be formed. If we instead assume that a grain surface 
containing no atoms of X and more than two atoms of Y 
will always be a site of Y 2 formation, Eq. (|23|l should be 
replaced with 

Kv Y + OL Ct , 

r(Y a « i o 9 P0 '° 28 

In the case of a = 0.3, 5 = 0.1, r = 10, this results 
in an improvement of a factor of about 2.4 in the value 
obtained for the production rate of Y 2 . This results in X 2 , 
XY, and Y 2 abundances, calculated with Eqs. H25fl to l|27|) . 
which differ from the exact results by 42%, 3%, and 26%, 
respectively; for X 2 and Y 2 this is a vast improvement 
from the results obtained when T(Y 2 ) was calculated with 
Eq. l]2.'{|l. The inaccuracy at high values of r also causes 
an inaccuracy in x(Y 2 ) at high (iV(Y)), which can be seen 
in Fig. |SJ, and probably contributes to the inaccuracies 
seen in x(X 2 ) and a;(XY) in Figs. and[SJi, at least at 
higher values of (N(Y)). In Fig.Efe-i, one can see that the 
results given by this approach for x(X 2 ) are inaccurate for 
low (7V(X)), at all values of < N(Y)>; at higher values 
of (N(Y)) the results are inaccurate for a larger range of 
(iV(X)). The results for x(XY) are also inaccurate for low 
(N(X)) for a range of (N(Y)). 

Fig-Elshows that the inaccuracy in the results given by 
the master equation method with low truncation values is 
high for high a. In particular, the inaccuracy was found 
to be high in the case of high a, low 5, and low r. This 
corresponds to the case of a high rate of accretion, a low 
rate of hydrogen desorption, and high rates for Reactions 
(j2Jl and J2J. This is the case in which there is a large 
population of X atoms on the surfaces of the grains. The 
inaccuracy is also high for high a, high <5, and high r, 
which corresponds to fast rates of accretion and hydrogen 
desorption, and a slow rate for Reaction J2Jl. This results 
in a large population of Y atoms on the surfaces of the 



grains. The results given by the master equation method 
with low truncation values are therefore inaccurate when 
there are large numbers of atoms on the surfaces of grains, 
as should be expected; this can also be seen in Figs. EJ-1 
andEJ-1. 

Fig. H shows the regions of (N(X))-{N(Y)) space in 
which the standard rate equation approach, or the low- 
truncation master equation approach, is more accurate. 
Plot (a) is for X 2 , plot (b) for XY, and plot (c) for Y 2 . 
Dark regions indicate that the rate equation approach 
is more accurate; light regions are those in which the 
low-truncation master equation approach is more reliable. 
Regions of intermediate shading indicate that both ap- 
proaches give results which are accurate to within 10%. 
As expected, when either the results given by one method 
or both methods are not accurate to within 10%, the low- 
truncation master equation approach is more accurate at 
low values of (iV(X)) and (N(Y)), and the rate equation 
approch at higher values. The better of the two approaches 
gives results for cc(X 2 ), x(XY), and ir(Y 2 ) which are usu- 
ally within 10 or 20% - and always within 25% - of those 
given by the exact method, except for the results given for 
.t(Y 2 ) in a small region of a—S—r space around a = 0.7, 
0.1 < S < 0.3, r > 0.7. In these cases, the more accurate 
method is the rate equation approach, and the value of 
x(Y 2 ) given by that method is around 40% less than the 
exact value. 

Although in general the rate equation approach is 
more accurate than the low-truncation master equation 
approach for higher values of (N(X)) and (N(Y)), it can 
be seen in Fig.0that for (N(X)) ~ 1.2 - 1.4, (N(Y)) ~ 3, 
the low-truncation master equation approach is more ac- 
curate for X 2 . However, in most of these cases the results 
given by the two methods are very similar, and the dis- 
crepancies between them and the exact results are within 
about 20%. 

4. Conclusions 

We have studied grain surface chemistry with a variety of 
approaches with a range of parameters. We deliberatelty 
made our study as general as was reasonable given the 
amount of work involved in a multiparameter study. We 
did so in anticipation of revisions of rates pertaining to 
specific problems treated in earlier literature, and of future 
studies of surface astrochemistry in a very wide range of 
environments including some in which atomic hydrogen 
accretion is unimportant. 

The modified rate equation method was found to be no 
more accurate than the standard rate equation approach 
under the conditions examined; indeed it was in many 
cases less accurate. The modified rate equation method 
that we used is closely related to that of Caselli et al. 
(2002). They studied a system in which reaction barriers 
are important, and the values of a approriate for their 
studies are large. The modified rate equation approach 
that we used gives reliable results for large a for some val- 
ues of r. Hence, our findings are in harmony with those of 
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Fig. 7. Regions where results for (a) X 2 , (b) XY, and 
(c) Y2, given by standard rate equation approach (darker 
regions) and low-truncation master equation approach 
(lighter regions) are more accurate. Regions of intermedi- 
ate shading are where both approaches give results which 
are accurate to within 10%. 



the specific a <C 0.1 and large 5 cases they examined. 
They, like the authors of a number of papers referenced 
here, were examining a region of parameter space they 
thought to be relevant when hydrogen accretion is impor- 
tant for the surface chemistry and the surface reactions 
proceed without barriers. Of the approximate methods 
we have examined, only the low-order truncated master 
equation approach is reliable for general small a, large 5 
cases. 

The approximate method of Green et al. (2001), which 
was devised specifically for the H,0 system, was also found 
to give inaccurate results under certain conditions, but can 
be improved by considering the possibility that an atom 
of species Y may accrete onto the surface of a grain on 
which there are already two Y atoms. 

The standard rate equation approach was found to 
work well except under conditions which lead to there be- 
ing a small average number of X and Y atoms on the 
surface of a grain. Under these conditions, the master 
equation approach with low truncation values is accurate. 
Therefore, a combination of the two approaches should be 
enough for a reasonably accurate calculation of the X2, 
XY, and Y 2 abundances in the parameter regime we con- 
sidered - generally within about 20% of the exact values. 
Practical approaches to the use of stochastic models are 
under investigation (Stantcheva et al. 2001). 
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Fig. 5. Percentage discrepancies in results given by the various approaches, plotted against (iV(X)} and (N(Y)) for 
all r. a - c: standard rate equation approach; d - f: modified rate equation approach; g - i: approximate method; j - 
1: low-truncation master equation approach, a, d, g, and j: discrepancy in x(X2); b, e, h, and k: discrepancy in x(XY); 
c, f, i, and 1: discrepancy in x(Y 2 ). Darker regions indicate greater discrepancies. 
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Fig. 6. Percentage discrepancies in results given by the various approaches, plotted against (iV(X)} and (N(Y)) for 
t = 1. a - c: standard rate equation approach; d - f: modified rate equation approach; g - i: approximate method; j - 
1: low-truncation master equation approach, a, d, g, and j: discrepancy in x(X2); b, e, h, and k: discrepancy in x(XY); 
c, f, i, and 1: discrepancy in x(Y 2 ). Darker regions indicate greater discrepancies. 



